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The computationally-determined limiting performance of shock 
isolation systems has been a useful tool in providing 
characteristics of optimal shock isolation. The limiting 
performance is defined as the minimum peak value of certain 
responses while other system responses are constrained. As 
is the case with most optimization problems, the 
"trajectory" in reaching the minimum performance index (peak 
response values) is unique, as is the minimum performance 
index itself. However, the responses of the system after 
the minimum performance index is achieved are not 
single-valued. This paper shows how unique isolator forces 
and corresponding responses can be chosen by super imposing a 
minimum settling time onto the limiting performance of the 
shock isolation system. Basically, this means that the 
system which has reached the peak value of the performance 
index is "settled" to rest in minimum time. 


INTRODUCTION 

The limiting performance of a system is its absolute optimal response 
characteristics. It is computed by replacing those portions of a system being 
designed by active generic isolator forces. These isolator forces are then 
obtained so as to minimize a given performance index while typically satisfying 
bounding constraints on response variables or isolator force magnitudes. Since 
isolator forces are not restricted to represent any particular design elements 
during the optimization procedure, the resulting limiting performance response is 
optimal over all possible design configurations. No conditions are placed on the 
number or type of elements which are replaced by isolators; they may be active, 
pasive, or nonlinear. For the class of problems treated in this paper, the 
performance index and the constraints are linear combinations of system response 
variables and isolator forces. Also, the equations of motion are linear, so that 
it is possible to formulate the optimization procedure as a linear programming 
problem. 

The limiting performance may be illustrated graphically by plotting the 
performance versus a constraint bound. If the performance index is chosen to 
minimize the maxmum response of a system subject to a prescribed constraint, the 
resulting tradeoff curve depicted in Fig. 1 gives the limiting performance of the 
system. The limiting performance charateristics are of considerable value to the 
mechanical system designer. First of all, they indicate from the design 
specifications alone whether a proposed design is feasible. Second, during the 
design cycle, they provide a measure of the success of the design configurations 
under consideration. Reference [1] describes limiting performance as applied shock 
isolation systems. Steady state systems [2], techniques for using general purpose 
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structural analysis computer programs to generate the equations of motion for 
limiting performance studies [3], and the use of limiting performance 
characteristics in identifying the optimum design of suspension systems for 
rotating shafts [4] have been treated. 



Fig. 1 Limiting Performance Characteristics 


Although the limiting performance provides useful information, it has been 
noticed that the min-max norm of the limiting performance gives a unique solution 
only until the peak value of the performance index is achieved. A non-unique 
solution occurs after the peak value. For the tradeoff studies, the response after 
the peak is of little importance as long as the unique performance index is 
obtained. However, rapid settling of the disturbed system due to the external 
disturbance is often desired. Therefore, it is necessary to impose an additional 
measure of performance to obtain a unique solution after the peak value of the 
performance index is achieved. The response after the peak is selected to achieve 
the minimum settling time. Two different approaches can be used to achieve this 
goal [5]. This paper deals with the formulation using the performance index and 
its application to the shock isolation problems. 


PROBLEM STATEMENT 

A linear vibrating system with n degrees of freedom subject to arbitrary 
external excitations f(t) and isolator forces u(t) is expressed in the first order 
system of differential equations 

s(t) = As(t) + Bu( t) + Cf(t) (1) 

where s(t) is an n-dimensional state vector and A, B. and C are n x n, n x nu and n 
x nf constant coefficient matrices. The quantities nu and nf are the number of 
isolator forces and excitations, respectively . Constraints are imposed on the 
dynamic system under study. The format of the constraints is 

X L < QjS + QgU + C^f < Xy for t Q < t < t f (2) 

where y^ and y^ are nc-dimensional lower and upper constraint vectors; Qj , Qg. and 
Q are nc x n, nc x nu, and nc x nf constant coefficient matrices; and t and t f 
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are the given initial and final times. 

The problem is to find an optimal isolator force u(t) which will transfer an 

initial state s(t ) = s to a desired final state s(t r ) = s~ l n the minimum time 
o — o i I 

while extremizing a given performance index of the form 

maX qn *-p rp 

Minimize J = {t Q < t < t^, IejS + EgE + E 3 — I) (3) 

where Ei • Eo* Eo are given n - nu > suid constant coefficient vectors. Since 

i. «3 

the min-max norm of the limiting performance gives a unique solution only until the 
peak value of the performance index is achieved, an additional measure of 
performance is desired to obtain a unique solution after the peak value. The 
resulting unique solution is referred to as limiting- performance/minimum - time 
(LP/MT) solution. 


LINEAR PROGRAMMING FORMULATION 

To obtain the LP/MT solution, the performance index given in Eq. (3) is 
modified. Two sets of performance indices are considered. One set of them, 
referred to as the transient performance index, is given by 

ItdX m fp 

J t = t Q < t < t t |eJs + Eg!! + E 3 il ( 4 ) 

where t is the time limit for the transient period. The other set, referred to as 
the steady-state performance index is defined as 

J s = t t < t < t f |e[s + Eg!! + Egil ( 5 ) 

Now, the ’’global" performance index is defined by 

j = j t ♦ -> s <6) 

Note that the vectors gj, E 2 * 2111(1 % are not changed in Eqs. (4) and (5). 

To place the optimization procedure into the standard linear programming form, 
the system in Eq. (1) is discretized using uniform time intervals to obtain a set 
of state difference equations 

s(k+l) = Gs(k) + H[Bu(k) + Cf(k)] (7) 

where s(k) = state vector at time t = t^ 

u(k), £(k) = isolator force and external 
assumed to be constant over 

G = e* h 

H = / h e* (hT) dr 

J o 


excitation vector at t = t^, 
the Interval t. < t < t. , . 
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h = time step = t^ + j - t^ (k=l , 2 N-l) 

The state vector, at any time t = t^, can be expressed as a function of the initial 

state s(l) and the isolator force history u(l), u(2),..., u(N-l) and the external 
excitation f(l) , f(2), .... f(N-l). For k = 1,2, ... ,N-1 

, k-i . . 

s(k+l) = G K s ( 1 ) + 2 G k J H[Bu(j) + Cf(j)] + H[Bu(k) + Cf(k)] (8) 

j=l 


The constraints in Eq. (2) are discretized similarly 

Y L (k) < Q x s(k) + Qgufk) + Q 3 f(k) < ^(k) for k = 1.2 N-l (9) 

The objective functions of Eqs. (4) and (5), which reflect the min-max norm, are 
discretized and converted into a constraint set. Since J^. is the maximum value of 

T T T 

|p!s + piu + p~f I for t < t < t. and so is J for t. < t < t., 

1“ Z" o O — t St I 

|pjs(k) + pju(k) + P 3 f(k)| < J t for t Q 1 t < t { 

IPjS(k) + p^uCk) + p 3 £ (k) | < J s for t t < t < t f 

To place this optimization problem into a standard linear programming form, 
def ine 


z 



( 11 ) 


where u = [ u(l)^ u(2)^ ... u(N-l)^ (12) 

and 

c T = [ 1 1 0 0 ] (13) 

Then the linear programming problem is to minimize 

J = c T z (14) 

subject to the constraints 

Hz < b (15) 

where H and b represent constraints of Eqs. (9) and (10). 

The minimum time (t . ) is the smallest time which will make the global 

v min J 

performance index of Eq. (6) stay within a desired value. Since the performance 
index can be computed for each iteration, an interpolation method such as the 
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secant method or simple bisection method [6] can be employed to find t 
efficiently. 


NUMERICAL EXAMPLES 

Example 1: A Single DOF System Subject to a Shock Velwave 

A single degree-of-f reedom (DOF) system composed of a mass m and supporting 
structure (Fig. 1) is subject to the horizontal shock velwave of Fig. 2. 



Fig. 1 A Single DOF Syatem 



Note: Slope (S), Time of Arrival (TOA), and Total =1+2 
Fig. 2 Shock Velwave 


Suppose the acceleration of the mass is to be limited to 15g (g = acceleration of 

x 

gravity). The optimal isolator force u (t) which minimizes the rattlespace between 
the mass and the supporting structure is desired. 

The equation of motion is 

z = u/m = U (16) 

The system is assumed to be at rest initially. The performance index to be 
minimized is 

max |z - y| (17) 
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where 


y 



v(t) dt 


and the constraint on the peak acceleration would be 


(18) 


|z| < 15g 

(19) 

Define a state vector 


. -p 

s = [ z z y ] 

(20) 


The problem can now be transformed into the standard LP/MT isolator problem format 
described previously. The optimal isolator is sought, which will reduce the 
disturbed rattlespace to zero in minimum time while minimizing the performance 
index and satisfying the constraints. The resulting time responses are shown in 
Fig. 4. The performance index is 0.914 in and the minimum time is 0.17 sec. 

Example 2: Two DOF Model of a Flexible Package Structure 

A two mass model of a flexible package structure with a rigid base is shown in 
Fig. 5. 


x 2 



Fig. 5 Two DOF Flexible Package Model 

The base is subject to external displacement which is described by 

f(t) = 12t 2 e _t [in] (21) 

The optimal isolator force is sought which will reduce the absolute displacement of 
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The performance index is to minimize 


T max | » | 

J = t l z 2 I 

while satisfying prescribed constraints 

|u| - lu/mjl 1 

1*2 1 - K -*l I **2™ 

Kl - K ' f l 4 X lna* 

The equations of motion can be written as 


(22) 


(23) 


]17. X + X(z x - z 2 ) = U ( 2 4) 

z 2 + \(z 1 - z 2 ) = 0 

where 

v = m 1 /m 2 (25) 

X = k/m 2 

Let 

S = [ Zj Zj z 2 z 2 ] T (26) 

Then the equations of motion, the performance index, and the constraints as given 
by Eqs. (24), (22), and (23), respectively, can be converted into the standard 

LP/MT format. Choose U [nay = 2g [in/sec 2 ], X lmgx = 2 [in], X 2max = 1 [in], ft = 

o 

100, and X = 6.28 [rad/sec ]. The solution for the optimal LP/MT isolator shows 
that the performance indices are J = 5.729 [in] and J = 0.099 [in], and the 

t S 

minimum time is t . = 12.0 [sec]. Figure 6 shows the time responses. 


CONCLUSIONS 

The objective of this study was to show how unique isolator forces and 
corresponding responses could be chosen by superimposing a minimum settling time 
onto the limiting performance of the shock isolation systems. The limiting 
performance / minimum time characteristics were computed by linear programming. It 
was demonstrated that the super imposition of minimum settling time provided not 
only the value of the optimal performance index but also the minimum settling time 
which, in turn, gives unique solutions for shock isolation problems. The optimal 
LP/MT isolator characteristics can be used to check the feasibility of proposed 
design requirements and to measure the success of a given design during the design 
process by comparing the response of the designed system with that of LP/MT 
characteristics. Furthermore, the LP/MT charateristics would provide with a 
designer an insight to build a near optimal shock isolation system. 
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